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Abstract 

We investigate by numerical simulations the behavior of the power dissipated 

in a resistive load capacitively coupled to a Josephson flux flow oscillator and 

compare the results to those obtained for a d.c. coupled purely resistive load. 

Assuming realistic values for the parameters R and C, both in the high- and 

in the low- Tc case the power is large enough to allow the operation of such a 

device in applications. 
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I. INTRODUCTION 



Low noise measurements in the sub-millimeter range, e.g.^ in radioastronomy, require 
a stable, low noise local oscillator for the receiver. The Flux Flow Oscillator (FFO), a 
device made with a long Josephson junction, has been proposed as a good candidate for 
this application [0]. Some of the most important features of such a device are the following 
the output power is relatively large (0.1 — l^iW), the oscillator can be easily tuned 
in a wide band (75 — 500 GHz) by varying the bias and a magnetic field, and the emitted 
signal has a very narrow linewidth (130 kHz at 70 GHz [Q, less than 2.1 MHz in the band 
280 — 330 GHz P). As the signal generated by a local oscillator has to be coupled to a mixer 
or to a transmission line, in the literature different couplings to a load have been realized and 
studied P,|5HlO|. In some works, in particular, the FFO has been d.c. coupled to a small 
junction acting as a detector, or, possibly, as a mixer; in this case, the real drawback is 
that the junctions cannot be biased independently. Capacitive coupling not only overcomes 
this drawback, but it also allows to increase the power transferred to the load, since it 
eliminates the d.c. loading of the oscillator. In this paper we investigate and compare these 
two coupling techniques for an FFO. We note in passing that ours is not the first work to 
propose the use of a d.c. block between oscillator and load device-see, e. 5'., However, in 
contrast with |^, which described a complete and detailed integrated subsystem, we focus 
specifically on the effects of the coupling element on system performance. 



II. THE PHYSICAL DEVICE - THE FLUX FLOW OSCILLATOR 

The FFO consists of a long Josephson junction biased by a d.c. current Ib and driven by 
the effect of a magnetic field H^,, perpendicular to the length L of the junction and parallel 
to the barrier, into a dynamical state in which the unidirectional motion of flux quanta 



(Josephson current vortices, or fluxons) takes place . With reference to the configuration 
sketched in Fig. |l], the fluxons continuously penetrate from the left edge of the junction and 
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propagate to the right; this regime corresponds to the appearance of typical branches in 
the I — V characteristic. An array of vortices travels in the junction with phase velocity 
u = {Vdc/cdfiQHe)c, where Vdc is the average voltage across the junction, He is the external 
magnetic field applied in the —y direction, c is the propagation velocity of electromagnetic 
waves in the junction and d is the effective magnetic thickness of the barrier. The largest 
signal in the flux flow regime is obtained biasing the junction at the top of the flux flow branch 
in the I — V characteristic; the corresponding voltage is determined by the velocity matching 
condition: the velocity u approaches the Swihart velocity c and, hence, Vdc — (cdfio) H^. 
One should take into account that a possible self magnetic field and the focusing effect of 
the external magnetic field could make Vdc larger. We shall assume that in the stationary 



state the expression for the junction voltage is given by the Fourier series |12 



v[x. 



t) = Vdc + Y^ Vn{x) exp{jnuFFt) , (1) 



where the coefficients V^ are complex amplitudes and the Yl'^ denotes summation from 
n = — oo to n = +00, but without the term n = 0, and upp = (27r/$o) V^c is the fundamental 
angular frequency of the series. Since the amplitude of the fundamental component of the 
series (|Tp is dominant, the frequency of the FFO signal is 

fpF = Vdc/% . (2) 

This is just the frequency expected from the Josephson relations when a d.c. voltage Vdc is 
applied across a junction. 

An important parameter for an oscillator is the frequency range in which it can generate 
a signal useful for applications, or, in other terms, its effective bandwidth. Since, under 
the velocity matching condition {u — * c), the frequency of the FFO is proportional to the 
magnetic field, fpp = (cdfio/^o)Hc, the effective bandwidth of the FFO is bounded to the 
values of the magnetic field within which the flux flow dynamic state is stable. Its lower 
limit is given by the minimum value of magnetic field needed to have penetration of fluxons 
into the junction, i.e., by the critical field He mm = "^jc^j- The upper limit can be estimated 



by supposing the maximum voltage in the flux flow regime to be of the order of the voltage 
gap (or of the order of -R„/c in the high-T^ context), as the nonlinear internal dynamics is 
strongly attenuated above that value. The tuning of the FFO to the desired frequency is 
achieved by changing the magnetic field to move the flux flow branch and by setting the d.c. 
current to bias the junction near the top of the branch. 

III. MODEL AND COMPUTATIONAL TECHNIQUES 

In Fig. |T] there is a schematic representation of a long Josephson junction, the FFO, that 
produces a signal coupled to the load through a capacitance. For this investigation the load 
device is assumed to be purely resistive. This assumption is made with the hypothesis that 
such a load model is sufficiently accurate to give information on the practical advantages 
of the capacitive coupling in comparison with simple d.c. coupling; for other loads, an 
appropriate model could be chosen and analyzed using the same approach. 

The mathematical model describing the flux-flow oscillator is the perturbed sine-Gordon 
equation (PSGE), which in normalized form is 

•^xx - '■Ptt - simp = aift - P^Pxxt - 7 > (3) 
with the boundary conditions 

<f,{0,t) + (3ip,t{0,t) = -7] , (4a) 

ip^{L,t) + (3ip^t{L,t) = -T] - iL{t) . (4b) 

Here, ip is the phase difference between the junction electrodes, x is the spatial coordinate 
normalized to the Josephson penetration length Aj, t is the time normalized to the inverse of 
the Josephson angular plasma frequency uj j, the subscripts indicate partial differentiation, 
the term in a represents shunt loss due to quasiparticle tunneling (here assumed ohmic), the 
term in /3 represents surface loss in the junction electrodes, 7 is the distributed bias current 
Jb normalized to the critical current density jc, L is the normalized junction length, rj is the 
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external magnetic field in the plane of junction and perpendicular to its long dimension in 
units normalized to jc^j- In the R — C load network is the current normalized to jc^jw 
(where w is the width of the junction), Rl is the load resistance normalized to the (linear) 
characteristic impedance of the junction Zq, Cl is the coupling capacitance normalized to 
the capacitance Cq = I/ujZq and ui = 1/RiCl is the load angular frequency normalized 
to ojj] accordingly,the equation for the current i/, is 



IV. NUMERICAL RESULTS 

Eqs. (^ 1^) with the boundary conditions Eqs. (^, |4BD have been integrated numerically 
using a 4th order Runge-Kutta algorithm on a spatially discretized counterpart, varying R^ 
from 0.1 to 10^ and Cl from 10^^ to 10. In the following we list the most relevant results. 
We note in passing that Fig. |1| depicts schematically a sandwich-type tunnel junction struc- 
ture. Although the fabrication technology for high-Tc tunnel junctions is not yet completely 
mature, significant progress in this direction is being made [|l^]. Accordingly, in what fol- 



lows we have chosen model parameter values that are aimed in the direction of describing a 
lightly-hysteretic tunnel junction. 

In Fig. ^ we report the Fourier spectrum of the voltage at the right edge of the junction 
v{L,t) calculated over more than two hundred periods. The spectrum consists of a d.c. 
component, whose height is the average voltage Vdc a fundamental line at the frequency fpF, 
having amplitude Kic(/ff), and a number of harmonics, up to the 7th; the small background 
is due to the finite sampling, to the finite integration time and to numerical noise. We remark 
that the fundamental frequency turns out to be just equal to what one would calculate from 
the height of the d.c. component (taking into account the normalization), in full agreement 
with Eq. (0). Moreover, one sees that the fundamental harmonic is dominant and the 
amplitude of the other harmonics decreases exponentially with increasing order, as can also 



5 



be inferred qualitatively from the inset in the same figure. We did not calculate the voltage 
spectrum for every value of the load used in our simulations; rather, we performed for most 
of them a simpler test to check that the series could be truncated to the second term. 
In fact, we calculated the rms amplitude of the fundamental harmonic VacK^fpF)! from 
the signal f (L, t) using the formula for the Fourier series coefficient and the rms value of the 
a.c. component of the signal Vac{L,t) = v{L,t) — V; the comparison shows that they are 
always equal within 1%, so that in practice all the power is in the fundamental harmonic. In 
other words, in the case of the R — C load the flux flow signal Vac{L, t) maintains, with good 
approximation, a sinusoidal form for every value examined of Rl and Cl- Therefore, the 
output power Pl is assumed to be the power of the flrst harmonic of the signal transferred 
from the flux flow oscillator to the resistive load Rl- For our R — C loading network it is 
given by 



Pl 



1 + {uL/uJFFf 



(6) 



2Rl 

Here Pl is normalized to the Josephson power P^ = VJ/Zq, where Vj = ^qujj/2tt is the 
normalizing Josephson voltage. 

We compared the output power for the R — C load to that of d.c. coupling to a resistive 
load, which has been already studied by Zhang [Q, to emphasize the different behavior. 



We have, in our model, the pure-R load case by setting ul = in Eq. (||) and Eq. (||). For 
the sake of consistency, we used the same parameters of Ref. [Q, i.e., a = 0.25, /? = 0.005, 



L = 20, biasing the junction with 7 = 1.25 on a flux flow step obtained with 77 = 4; the 



results essentially agree with those in [|TJ], differences being attributable to the fact that 
with purely resistive coupling the output waveform is no longer a clean sinusoid, so the total 
power differs from the flrst-harmonic power. The maximum output power dissipated by the 
a.c. component is PLmax = 0.76 for Rl = 7, and changes very little in the range 3-10 (we 
note parenthetically that these numbers are reminiscent of those obtained in an early study 
of a resonant-fluxon oscillator with pure-i? loading). This value is represented in Fig. |^ 
and Fig. ^ as a horizontal straight line for the sake of comparison with the R — C case. 
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The R~C case is, at first sight, more complex, as has already been noted by Zhang (see 
Section 4.4 of Jl^)- Iii Fig- 3 we report the output power as a function of the resistance Rl for 
fixed values of the ratio ojff/^l- We see that the power increases as the oscillator frequency 
grows with respect to the load characteristic frequency, and that it is not a monotonic 
function of the resistance. In fact, maximum output power is obtained for ~ 1, i.e., when 
the load resistance Rl is close to the junction characteristic impedance Zq. This result is 
consistent with what one should expect from considering the best matching condition in the 
framework of (linear) microwave transmission line theory. We remark also that the variation 
of the output power is negligible for ujff/^l ^ 5. In Fig. ^ the comparison with the curve 
with Pf = 0.76 shows in which range the R — C loaded FFO is more efficient than the R 
loaded FFO. 

In Fig. ^, in order to provide a more straightforward tool for device design, we plot the 
output power as a function of the load capacitance for fixed values of the load resistance 
Rl. For a given value of Rl-, the power is first an increasing function of Cl which, however, 
quickly saturates. The asymptotic values lie on the curve with uff/^l = 100 of Fig. |^, 
which well approximates the case of infinite capacitance. 

Whereas, as mentioned above, the R—C case appears at first to be more complicated than 
the pure-i? case, in fact it is simpler: once the values of Vac and uoff '^'^ Eq- i| are established 
by the junction dynamics, at least for the parameter values that we have studied, the 
oscillator behaves with respect to the load much as a simple linear oscillator, characterized 
by a fixed internal voltage and a fixed internal resistance equal to Z^. This operational 
simplicity is obtained essentially from the elimination of the d.c. loading of the oscillator, 
which, instead, is present in the pure-i? case. 

V. DISCUSSION 

As we have seen in the previous section, in order to maximize the output power for a 
given value of the resistance, it is recommendable to increase the load capacitance, because 
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the best coupling is obtained for > 1. To check whether this indication can be translated 
into realistic parameters to design and fabricate a device, we shall estimate what one should 
expect in two practical cases. In both cases we shall estimate Cp from the formula Cp — 
€reoS/d, with the usual meaning of the parameters. First, we shall consider the high-Tc case, 
and we shall suppose that the coupling capacitor is a typical silicon monoxide element (e^ = 
5, d = 600 nm and S = 50/xm x50/im). Since cuj — ^j2T:jc/ {^qC) and Zq = {1/w)^ ii^d/C 
where C is the junction capacitance per unit area, the previous definitions give: 



Assuming from the literature jc — 20 kA/cm^, w — 200 nm, d — 280 nm, C = 30 fF///m^, 
we finally get Cl = 15.7. To evaluate the power dissipated in the load, we can use the same 
figures, and find that Pu — 0.13/xW; this gives an idea of the order of magnitude of the 
power that could be extracted. 

Next, let us consider a low-T^ oscillator. In this case, the coupling capacitor can be made 
by growing a niobium oxide film by anodization (e^ = 25, o? = 100 nm). From the previous 
formula, assuming typical parameters for Nb-A103;-Nb Josephson junctions {jc — 1 kA/cm^, 
w = 3 /im, d = 80 nm, C = 60 fF /fim'^), we find Cl = 2.3. Of course, in this case one 
is left with the problem of fabricating a load resistor that should be perhaps 0.1 ohm (or 
less), but this is feasible with present-day technology; moreover, our numerical simulations 
should be considered as being merely indicative, in that the dissipation in a low- Tc junction 
can be significantly lower than the value considered in this paper, which might bring about 
qualitative changes in the dynamics. Nevertheless, we think that it is interesting to take 
into account a design based on the well established niobium technology, and we shall explore 
this topic in more detail in the future. 
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FIGURES 



FIG. 1. Schematic layout of fiux-fiow oscillator capacitively coupled to a load. 

FIG. 2. Fourier spectrum of the output voltage. The waveform is shown in the inset. The 
parameters used in this calculation are: Rl = 1.5, Cl = 0.01, V^c = 4.218, rj = A, a = 0.25, 
/? = 0.005, 7 = 1.25, L = 20. The amplitude is normalized to Vj = ^QiOj/2'K; the frequency is 
normalized to fpF = 0.671. 

FIG. 3. Output power as a function of load resistance for different values of the ratio between 
the FFO frequency and the R — C characteristic frequency. Here r) = 4, a = 0.25, /3 = 0.005, 
7 = 1.25, L = 20. 

FIG. 4. Output power as a function of load capacitance for different values of the load resis- 
tance. Parameters are as in Fig. 3. 
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